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Introduction 



GALA (jMucciarelli et al.ll2013l ) is a code written in standard Fortran 77 and aimed at finding 
the best atmospheric parameters and the abundance of individual elements by using the equivalent 
widths (EWs) of metallic lines, providing graphical and statistical tools to evaluate the goodness 
of the solution. The derivation of the abundances is performed by using a mo dified version of 
the WIDTH9 cod^ (originally developed by R. L. Kurucz) in its Linux version ( Sbordone et al 



2004') ■ In the current releas e, GALA can manage the classical grid of ATLAS9 models computed 



(Meszaros et al 



12012 



bvlCastelh fc Kurucz (|2004l ). the grid of new ATLAS9 models computed for the APOGEE survey 



and the MARCS models with the standard composition (jGustafsson et al 



20081 ). When the ATLAS9 models are used, new model atmospheres are calculated starting from 
an existing guess model and according to the pre-tabulated Opacity Distribution Functions (ODF) 
and Rosseland opacity tables. When MARCS are used, each new model is obtained by interpolating 
within the MARCS grid. 



2. Installation 



Once you have downloaded the archive file gala. tar. gz from the website 



www. cosmic-lab. eu/Cosmic-Lab/Products .html 



type the commands 

gunzip gala. tar. gz 
tar -xvf gala. tar 

These commands unpack the archive file creating a directory named GALA/. Within this direc- 
tory there are three subdirectories ( src/, bin/ and tutorial/) together with some information 



^http://wwwuser. oat. ts.astro.it/castelli/sources/WIDTH. html. 
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files. The directory src/ includes all the source files and the Makefile needed to install GALA. The 
directory bin/ includes the executables created after the installation. In the directory tutorial/ 
there are some examples of the input files for GALA as reference. 

(1) Install ATLAS9 

Before starting the installation of GALA, you need to install in your machine the ATLAS9 code if 
you plan to use these models (obviously, if you plan to use only MARCS models you can avoid to 
install ATLAS9). GALA includes a dynamical call to the last version of this code, available on the 
website of F. Castelljl. Download the source code atlasQmem.f or with the command 

wget -r -nd http://wwwuser.oat.ts.astro.it/castelli/sources/atlas9/atlas9mem.for 

and compile it 

ifort -double-size 64 -save -o atlas9mem.exe atlas9mem. f or 

Note that GALA assumes that you have this code installed and that the executable is named 
atlas9mem.exe. Please do not compile ATLAS9 with different names. 

(2) Set the Makefile 

Firstly, you need to properly update the Makefile. GALA can be compiled with the Intel For- 
tran Compiler, that you can download from the Intel websitfl after the registration. Also, the 
graphical output files are created by the program galaplot using the SuperMongc0 (SM) libraries 
(namely libplotsub . a, libdevices . a and libutils . a, compiled in single precision) and the Xll 
libraries. Set in the first lines of the Makefile the path of your Intel Fortran Compiler and that of 
the directory where the SM libraries are located. 

(3) Set some useful paths 

Now, you need to properly set in the GALA/src/Paths . f file the paths of some directories used by 
the code: 

(1) 'exe' = it is the path where the executables of GALA are located (the subdirectory bin/); 

(2) 'at9exe' = it is the path of the directory where the ex ecutable atlas9mem . exe is located; 



(3) 'at9fc' = it is the path where the ATLAS9 models bv lCastelli Sz Kuruca (j2004l ) are saved, see 



http://wwwuser.oat. ts.astro.it/castelli/sources/atlas9codes. html 
^http: / /software. Intel. com / en- us / non-commercial-software-development 
^http://www.astro. princeton.edu/~rhl/sm/ 
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Section 12 

(4) 'at9apo' = it is the path where the ATLAS9 models by iMeszaros et alj (|2012l ) are saved, see 
Section 12 



(5) 'rn arcs' = it is the path where the MARCS models with the standard composition bv lGustafsson et al 



(|2008l ) are saved, see SectionlT) 

The directories of the ATLAS9 models need a precise structm'e. In each directory you need to 
create three subdirectories, namely ODF/, MODELS/ and MODELS-new/, the first including the ODFs 
and the Rosseland opacity tables, the second including the grid of model atmospheres used as guess 
models, while the third is the directory where the new models created during the run of GALA will 
be saved. The directory for the MARCS models does not need a specific substructure. Examples of 
these subdirectories, as well as more details about the use of the different model grids, are discussed 
in Sect. 13 

Note that you do not need to have necessarily all these models to use GALA. If you plan to use 
only a given kind of models, you can delete or leave blank the paths of the other models in the 
Paths. f file (but in this case do not try to call them when you run GALA). 

(4) And now...instan GALA! 

To install GALA, type the command 

make all 

and the executables will be saved in the subdirectory /bin. 

Finally, put the path of the subdirectory bin/ in your login file according to the shell environment 
of your machine (for instance in the configuration file .bashrc or .tcshrc). The directory /tutorial 
includes examples for the three input files. You can check the installation of GALA by using these 
files and simply type gala. 

Note that GALA is dimensioned to manage up to 1200 spectral lines for each star (parameter INL 
in the source file Declare. f). If a given star has more than this number of measured transitions, 
GALA skips this star providing a warning message. For our experience this value is reasonable 
but if you think the this limit is too low for your stars, you can increase the parameter INL and 
recompile the code (note that a uselessly high value for INL will lead to a waste of computer 
memory) . 



3. Installation problems 

We list some common installation problems (typically related to the combination Fortran 
compiler /system architecture) and their possible solutions. 

• The SM libraries compiled in double precision can give problems during the compilation or 
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(when the compilation weh works) galaplot creates corrupted or badly drawn postscript 
files. We recommend to use for galaplot SM libraries compiled in single precision. Because 
it is usual to compile SM in double precision, we suggest to perform two distinct installations 
of SM, in order to save wherever you want the libraries compiled in single precision. The 
instructions to perform the installation of SM are explained in the sm. install file in the 
directory of the SM source code. 

• A warning message appears during the compilation of galaplot like 

Id: warning in /Users/alessio/MAGAZZINO/SMONGO-lib-ELE/libplotsub. a, file is 
not of required architecture 

In this case, you can try to compile galaplot telling the compiler to generate the code 
for a specific architecture. You need only to add the option '-m32' (if your machine is a 32 
architecture) or '-m64' (if your machine is a 64 architecture) in the Makefile 

$ (FORTRAN) -m32 -traceback - galaplot galaplot. f mplotl.f -1X11 -L$(SMLIB) 
-Iplotsub -Idevices -lutils. 

• The SM commands under Fortran are unrecognized, with warning messages like 

Undefined symbols: "_sm_limits_" , referenced from: _mplot_ in if ortLFDlnn. o 

This happens because some incompatibilities between the different names of the SM com- 
mands as defined by Fortran and C compilers can occur. The users with experiences with 
other tools devoted to spectroscopic analysis (as MOOG and DAOSPEC) have already dealt 
with these problems. We refer the reader to the Chapter 4 of the DAOSPEC Cookbook as 
useful documentation to solve some of these problems. 

With the last versions of the Intel Fortran Compiler, the SM commands have an extra '_' as 
suffix (not needed if you compile galaplot with G77 Fortran Compiler for instance). The 
source file mplotl . f calls the SM commands as smJimits- and so on. Note that we provide a 
modified version of this file (named mplot2 . f ) where the SM commands are called with their 
standard names (i.e. smJimits); you can try to compile galaplot with this source file, if the 
error messages shown above appear. 

• If you cannot compile galaplot with the Intel Fortran Compiler because of incompatibilities 
between the SM libraries and the Fortran commands as called in the source code, you have 
two alternative options: (1) try to compile only galaplot with another Fortran compiler (as 



'''http://www.bo.astro.it/~pancmo/docs/daospec.pdf 
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G77 for instance), substituting the file mplotl.f witli mplot2.f; (2) if you cannot compile 
galaplot with other compilers, you can perform the installation of GALA without galaplot. 
In this case, you can easily disable the graphical option by commenting in the file gala.f the 
line 



call system ( 'galaplot ' ) 

and then compile the code with the command make gala. 



4. Executable files 



The executable files that will be created in the subdirectory bin/ are: 



1. gala — the executable file to run GALA; 

2. galaplot — this program produces the graphical output for each analysed star; 

3. interpoljnarcstoatlas — this is a modified version of the code to interpolate the MARCS 
model atmospheres developed by T. Masseron and available in its original version at the web- 
site http://marcs.astro.uu.se/software.php. The current version writes the output MARCS 
model in the ATLAS9 standard format; 

4. marcstoatlas — this routine reads an usual MARCS model atmosphere converting it in the 
ATLAS9 standard format. This routine is a modification of the program (developed by B. 
Edvardsson) available in the Uppsala website (http://marcs.astro.uu.se/software.php) and it 
can be used independently by GALA with the command-line instruction 



cat < model-name > / marcstoatlas 



(where < model-name > is the input MARCS model), storing in the Fortran unit 88 the 
ATLAS9-like output model; 

5. width9_gala — this is a modified version of the WIDTH9 code available on the website of F. 
Castelli|j. The main d ifferences are the formats of the input and output files. Some control 
cards of WIDTH9 (see ICastellilll988l ) have been disabled because not used by GALA; 

6. extract_kur — this program reads the classical linelists of Kurucz/Castelli stored in a file 
named master_linelist . tmp and a file (named input.dat) including your linelist (in the 



®http://wwwuser. oat. ts.astro.it/castelli/sources/WIDTH. html 
^littp://wwwuser. oat. ts.astro.it/castelli/linelists. html 
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first column the wavelength in A and in the second column the code of the element in the 
Kurucz formalism: i.e. 26.00 for Fel and 26.01 for Fell). The code matches the two files and 
writes in the output file (output.dat) the information needed to the GALA input. 

7. write_autof 1 — this program writes the basic structure of the input file autof 1 .param (de- 
scribed in Section [5. ip . with the proper name of each keyword. 



5. Input files 

• autof 1. param — it includes the main configuration parameters used by GALA to perform 
the analysis (described in Section 5.1); 

• list_star — it is the list of the stars to analyse and their guess atmospheric parameters 
(described in Section 5.2); 

• a file (one for each star listed in list_star) with extension .in, including the list of the 
measured lines, their EWs and their atomic parameters (described in Section 5.3). 



5.1. Layout of autof 1. param 

The input file autof 1. param includes two groups of keywords, namely the basic parameters 
and the subordinate parameters. The command write_autof 1 writes the basic structure of the file 
(with some parameters already filled by default values); note that the order of the keywords (as 
well as the presence or not of the header lines) does not really matter. A template of this file is 
available in the subdirectory tutorial/. 

(1) Basic parameters: these are the parameters that you need to set in order to perform the 
analysis: 

• model specifies the kind of model atmospheres; the accepted values are 

-atlas9 for the ATLAS9 model atmospheres available in the website by F. Castelli; 

-apogee.cmlSO AP0GEE-ATLAS9 grid with [C/Fe]=-1.50; 

-apogee_cml25 AP0GEE-ATLAS9 grid with [C/Fe]=-1.25; 

-apogee.cmlOO AP0GEE-ATLAS9 grid with [C/Fe]=-1.00; 

-apogee.cm075 AP0GEE-ATLAS9 grid with [C/Fe]=-0.75; 

-apogee.cmOSO APOGEE- ATLAS9 grid with [C/Fe]=-0.50; 

-apogee.cm025 APOGEE- ATLAS9 grid with [C/Fe]=-0.25; 

-apogee.cpOOO AP0GEE-ATLAS9 grid with [C/Fe]=+0.00; 

-apogee.cp025 AP0GEE-ATLAS9 grid with [C/Fe]=+0.25; 

-apogee_cp050 AP0GEE-ATLAS9 grid with [C/Fe]=+0.50; 
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# BASIC 


PARAMETERS 


mo dsl 


atlas 9 








g 


uncer^ 





qiron 


26. 00 


epot 





grc 





vtc 





metal 





stepgfr 


0. 1 


steptef f 


SO 


stepvt 


0. 1 


e™in 


-S.S 




-2.S 


error 


10 


iter 


10 


# SUBORDINATE PARAMETERS 


r] 


1 




3 


mimf e 


10 


ep Lever 


l.S 


maxsigf 


1. 


mini 


3 


errfit 


1 


par opt 





inters 





kuriter 


s 


plot 


1 


verbose 


1 



-apogee.cp075 AP0GEE-ATLAS9 grid with [C/Fe]=+0.75; 
-apogee.cplOO AP0GEE-ATLAS9 grid with [C/Fe]=+1.00; 

-marcs for the MARCS model atmospheres with the standard chemical composition. 

We refer to Section [7] to detail about the download of these models and the way to manage 
them in order to use correctly GALA. Alternatively, you can specify the name of a given 
model (in ATLAS9 format) |f| located in the working directory. If the specified model is found 
in the current directory, all the optimization options (together with the guess and refinement 
working-blocks) are switched off and only the analysis working-block is executed (with fixed 
parameters). Otherwise, if the specified model does not exist in the directory, GALA stops 
the procedure and exits with a warning message. 

• guess enables the guess working-block, before the analysis working-block. The parameter is 
the number of iterations performed during this working-block. Guess= (or negative values) 
disables the execution of this working-block. If you decide to use this option, we suggest to 
adopt at least 3 iterations for the guess working-block. 

• refine enables the refinement working-block, after the analysis working-block. The parameter 



*If you wish to use a model atmosphere calculated with the ATLAS12 code ( Castellill2005al ). you need to erase 
the 22 lines corresponding to the 'ABUNDANCE TABLE' control cards. Also, to use a specific MARCS model, you 
have to put it in ATLAS9 format with the program marcstoatlas. 
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accepts the values (disabling the working-block) and 1 (enabling the working-block). 

• uncer enables the calculation of the uncertainties in the derived chemical abundances due to 
the atmospheric parameters. The accepted values are 
: no uncertainties calculation is performed; 



1 : the uncertainties are calculated according to the prescriptions by ICayrel et al.l (j2004l ) , 
only if Tgfj has been optimized spectroscopically; 

2 : together with the uncertainties estimated with the previous method, also the uncertainties 
due to the atmospheric parameters following the classical method are calculated, varying 
each time one parameter only (keeping the others fixed) and without taking into account the 
covariance terms. 

qiron indicates the code of the element (with its ionization stage) whose lines are used to 
perform the optimization. The formalism is the same used in the Kurucz codes (i.e. 26.00 
for Fe I and 26.01 for Fe II). In principle, GALA accepts all the elements but obviously only 
a few elements (i.e. Ti or Ni) have a sufficient number of available lines to perform a reliable 
spectroscopic optimization. 

epot enables the spectroscopic optimization of the temperature. The accepted values are 

: the temperature is not optimized and rather fixed to the input value; 

1 : the temperature is found by erasing any trend in the A(Eljj — Xex plane (where Xex is the 
excitation potential in eV). 

grc enables the spectroscopic optimization of the gravity. 

: the gravity is fixed to the input value; 

1 : gravity is optimized in order to obtain the same abundance from neutral and single ionized 
iron lines; 

2 : if this option is chosen, GALA reads in list_star the term e = log(4GM7rcr/L) and for 
each iteration log g is computed as log g=e+4-Teff (see Section 5.2); 

3 : the gravity is computed with a second degree relation logg= A-|-B-Tefj-|-C-Tjj. In this 
case the coefficients A, B and C are provided in the list_star file (see Section 5.2). 

vtc enables the optimization of the microturbulent velocity. 

: the velocity is fixed to the input value; 

1 : the microturbulent velocity is found by erasing any trend in the A(E1) — EWR plane, 
where EWR is defined as log(£'Vl^/A), where EW and A are expressed in the same units. 

metal enables the optimization of the overall metallicity of the model atmosphere, that is 
chosen in order to reproduce the average abundance derived from lines of the element specified 
by the parameter qiron (within the metallicity step of the employed models grid). 



'A(EI) indicates for a given element El the number abundance as calculated by WIDTH9. 
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• stepgr the step used to investigate log g in the analysis and refinement working-blocks; 

• steptefF the step used to investigate T^s in the analysis and refinement working-blocks; 

• stepvt the step used to investigate Vfurb in the analysis and refinement working-blocks; 

• ewmin, ewmax these are the minimum and maximum value for EWR. All lines outside 
the EWR range defined by these parameters are excluded from the analysis. If you do not 
know a priori the optimal EWR range, a first run of GALA with a wide range including 
all lines is recommended (with EWR ranging from ewinin= -10 to ewniax= 0), followed by a 
more informed second run, with EWR limits deducted from visual inspection of the graphical 
outputs. 

• error the maximum error in percentage in EW, in the case in which the error for each line 
is provided in the linelist file. All the lines with error larger than this threshold are excluded 
from the analysis. 

• iter the maximum number of allowed iterations. This parameter avoids the risk of infinite 
loops. A reasonable value is 10 (if the input parameters are not so far from the real solution) , 
but you should adapt this value to her/his particular needs. 

(2) Subordinate parameters: these parameters specify some finer details of the analysis; we 
suggest to use the default values for these parameters, unless you have specific needs 



• rj specifies whether the line-rejection uses the mean (rj=0) or the median value (default 
rj= 1). 

• smax is the number of a used in the line-rejection process (default smax= 3). 

• minnfe corresponds to the minimum number of surviving lines that allows to perform the 

optimization. In principle, the spectroscopic optimization of the atmospheric parameters is a 
statistical method which requires a large (at least 40-50) set of lines. Optimizations derived 
with a low number of lines arc afi'ccted by fluctuations due to the small number statistics and 
the lines distribution in the A(E1) -Xex and A(E1) — EWR planes (but these uncertainties are 
taken into account in the calculation of the errors through the Jackknife bootstrap technique) . 
You should be careful that a very low number of lines (typically ten or less) would provide 
very dangerous results (default minnf e= 10). 

• eplever the spectroscopic derivation of T^e needs a large range of Xex covered by the used 
lines, in our opinion at least 1.5-2 eV. For smaller range, the optimization of Tgff can be 
unreliable. This keyword specifles the minimum Xex range (in eV) allowed to perform the 
optimization of T^fi (if requested) . If the range of Xex (calculated considering only the lines 



- 10 - 



surviving line-rejection) is smaller than eplever, then the optimization of Tcff is swichted off; 
eventually, if the other parameters are to be optimized, their optmization continues (default 
eplever= 1.5). 

• maxsig is the maximum dispersion by the mean or median for the element chosen with 
qiron allowed to run GALA. When the dispersion exceeds this value, all the optimizations 
are disabled and the analysis is performed with the input atmospheric parameters specified 
in the list_star file (default maxsig= 1.0). 

• mini corresponds to the minimum number of lines to perform the process of line-rejection. If 
the number of lines is less than the value of mini {Nunes <minl), the lines are rejected only 
on the basis of uncertainty in EWs {(Tew) 'ind of the pair of values ewmin-ewmax (default 
iiiinl= 3). 

• errfit specifies the kind of linear fit performed on the A(E1) — Xex and A(E1) — EWR planes. 

: the linear fits are computed without considering the uncertainties in x and/or y; 

1 : the linear fits are computed by taking into account the uncertainties in A(E1); 

2 : the linear fit in A(E1) — Xex is computed including the uncertainties in A(E1) and that 
in A(E1) — EWR including the unc ertainties both in EWR and A(E1), following the method 



implemented by lPress et al.l (|l992l ). 

The use of the uncertainties in the slopes computation is recommended, but only if you have 
reliable errors in EWs, able to provide a realistic ranking among the used lines. If you do not 
trust your EW uncertainties, use errf it= 0. Note that if errf it= 1 or 2 but at least one of 
the input lines has aEW= 0.0, the inclusion of the errors in the slope calculation is disabled 
(default errf it = 1). 

paropt specifies the kind of optimization parameter used for Teg and vt'- 

: optimizes these parameters by using the slopes (default); 

1 : by using the Spearman ranking coefficient. 

interz enables the interpolation to the zero value of the optimized parameters (interz= 1). 
This option is included to provide visually nicer results if one wants to plot temperatures and 
gravities in a T^g — logg diagram. It interpolates the slope in the A(E1)-Xea; and the A(E1)I- 
A(E1)II to zero, and determines the corresponding Tefj and logg, thus avoiding the grid effect 
in the results. If this interpolation is switched off, all results will be in steps specified by the 
stepgr, steptef f , and stepvt parameters, while if the interpolation is switched on, they will 
be distributed more continuously. Note, however, that all uncertainties will still be computed 
on the best slopes, not on the interpolated zero slopes (default interz= 0) . 

kuriter specifies the number of blocks of 15 iterations each used by ATLAS9 to calculate a 
new model atmosphere (default kuriter= 5). This option is ignored if the keyword model is 
referred to MARCS models. 
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• plot enables the graphic output in postscript format 

: does not produce a plot; 

1 : saves the plot of the results in a postscript file (default option) . 

• verbose specifies the verbosity level on the terminal 

0: no message at all. It advises only when GALA ends the analysis; 

1: the sequence of the analysed stars is shown, without additional information; 

2: this is like verbose= 1 but also specifies the main steps performed by GALA (without 

additional information about them); 

3: this is the default verbosity level. It is like verbose= 2 but also displays information about 
the line rejection, the derived abundances and concerning the Guess Working-Block. 



5.2. Starlist file 

The file list_star includes the sequence of the stars that you plan to analyse, together with 
the input atmospheric parameters. The structure of the file depends basically by the kind of the 
adopted optimization for the surface gravity (as specified by the basic parameter 'grc' in the file 
autof 1 .param). 

• if you adopted grc= or 1, each line of the file will be as follows: 



ngcl806_2E 4000 

ngclS06_27 4200 

ngcl806_29 4100 

ngclS06_30 4400 

ngcl806_32 4100 



0.7 -O.E ^. 

1.2 -0.5 2.0 

1.0 -O.E 2.0 

1.5 -0.5 2.0 

1.0 -O.E 2.0 















































































Fig. 1. — Example of the list_star file when the option grc=0 or 1 is selected. 



1. the first column contains the root name of the input files, without the extension . in; 



2. the four following columns contain the guess atmospheric parameters in the following 
order: T^s, log g, [M/H] and vt- Note that if you provide a wrong value of [M/H] (not 
available in the chosen grid of models) , GALA changes it and assumes the closest metal- 
licity; 

3. the sixth column contains the value of the a-enhancement adopted in the computation 
of the ATLAS9 models. Also in this case, if you provide a value of [ct/Fe] not available 
among the models, the closest value is used. If you use MARCS models GALA neglects 
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this parameter; 

4. the last three cohimns are the errors associated to Tcfj, logg and vt respectively. If 
the parameters are optimized spectroscopically, these uncertainties are ignored (you can 
also put for all of them, as in the shown case); on the other hand, if a given param- 
eter is fixed and the uncertainties with the classical method are requested (uiicer=2 
in autof 1 .param), these values are used to vary each parameter computing again the 
abundances. If in the input file the parameter uncertainties are but you request to 
calculate the abundance uncertainties, the steps used to investigate the parameter space 
(and defined by the keywords steptef f , stepgr and stepvt among the basic parameters 
in autof 1. param) are assumed as parameter errors. 

• if you adopted grc= 2, the input file will be 



ngclSOS 25 


4000 


0. 


7 


-0. 


5 


ngclSOe 27 


4200 


1. 


2 


-0. 


5 


ngclSOS 29 


4100 


1. 





-0. 


5 


ngclSOS 30 


4400 


1. 


S 


-0. 


E 


ngclSOS 32 


4100 


1. 





-0. 


5 


ngclSOS 40 


4100 


0. 


S 


-0. 


5 



2. 





0. 


G 











-13. 





2. 





0. 


G 











-13. 


2 


2. 





0. 


G 











-13. 


1 


2. 





0. 


G 











-13. 


2 


2. 





0. 


G 











-13. 


2 


2. 





0. 


G 











-13. 


4 



Fig. 2. — Example of the list_star file when the option grc=2 is selected. 



The additional tenth column e is calculated according to the Stefan-Boltzmann equation as 
e = log(4GM7r(7/L), where G is the gravitational constant, a the Boltzmann constant and M 
and L are the mass and the luminosity of the star. In this case the surface gravity will be 
calculated as logg= e + 4- log T^s 



• if you adopted grc= 3, the input file will be 



ngclS06_2S 4000 O.T -O.S 2.0 0.0 -T. 69087 0.00198 

ngcl80S_27 4200 1.2 -O.E 2.0 0.0 -T. 69087 0.00198 

ngcl806_29 4100 1.0 -O.S 2.0 0.0 -7.69087 0.00198 

ngcl306_30 4400 l.S -O.S 2.0 0.0 -7.69087 0.00198 

ngcl806_32 4100 1.0 -O.S 2.0 0.0 -7.69087 0.00198 



Fig. 3. — Example of the list_star file when the option grc=3 is selected. 



where A, B and C are the coefficients of the quadratic relation log g= A+B-Tes+C-T^f 



5.3. Linelist file 



For each star (specified by the <rootfilename> in the star_list file) an input file <rootf ilename> . in 
is needed. If the linelist file of a given star listed in list_star is not available in the directory. 
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GALA provides a warning message and moves to the next star. Fig. H] shows an example of a 
hnehst input file for GALA. Note that this file does not need a specific numeric format, but it must 
contain the following information for each measured line: 

• the wavelength (expressed in A); 

• the observed EW (expressed in mA); 

• the uncertainty in EW (expressed in mA); 

• the code of the element in the usual Kurucz notation; 

• the logarithm of the oscillator strength; 

• the excitation potential (expressed in eV); 

• the logarithm of the radiative damping constant, 7r-adE§; 

• the logarithm of the Stark damping constant, jstark] 

• the logarithm of the Van der Waals damping constant, ^vdw'^ 



the velocity parameter a as defined bv lBarklem et al.l (j200fl|) and usually available for many 
metallic lines both in the Kurucz/ Castelli linelist and VALD databases. If the parameter a 
is not available for a given line, the value can be substituted with 0. 



■5403 


0220 


124 


80 


4 


040 


26 


00 


-1 


03 


4 


076 


7 


84 


-5 


89 


-7 


73 





00 " 


5406 


7750 


65 


30 


3 


610 


26 


00 


-1 


46 


4 


371 


8 


24 


-4 


50 


-7 


19 





28 


S408 


9400 


216 


60 


11 


980 


22 


00 


-3 


84 





000 


6 


23 


-6 


16 


-7 


74 





25 


£410 


9100 


140 


70 


4 


££0 


26 


00 





40 


4 


473 


8 


10 


-£ 


29 


-7 


13 





28 


S412 


7860 


49 


80 


4 


040 


26 


00 


-1 


72 


4 


434 


8 


23 


-5 


16 


-7 


15 





28 


S42S 


2500 


89 


10 


5 


090 


22 


00 


-3 


01 





021 


G 


23 


-G 


16 


-7 


73 





25 


S429 


1380 


113 


30 


6 


080 


22 


00 


-0 


44 


2 


345 


8 


08 


-6 


04 


-7 


77 





00 


£436 


2960 


62 


10 


4 


230 


26 


00 


-1 


£1 


4 


386 


8 


24 


-4 


61 


-7 


19 





28 


S44S 


0420 


134 


20 


4 


220 


26 


00 


-0 


02 


4 


386 


8 


25 


-4 


86 


-7 


19 





28 


S4S3 


6430 


111 


40 


3 


540 


22 


00 


-1 


61 


1 


443 


7 


98 


-G 


10 


-7 


79 





00 


S460 


8730 


47 


30 


4 


410 


26 


00 


-3 


58 


3 


071 


7 


11 


-G 


20 


-7 


81 





26 


£461 


SSOO 


GO 


90 


3 


790 


26 


00 


-1 


64 


4 


446 


8 


27 


-4 


94 


-7 


16 





20 


£462 


4850 


SI 


00 


4 


110 


28 


00 


-0 


93 


3 


847 


8 


21 


-4 


06 


-7 


22 





28 


5463 


2760 


130 


80 


4 


930 


26 


00 





07 


4 


434 


8 


21 


-4 


97 


-7 


17 





28 


5473 


9010 


107 


50 


21 


710 


26 


00 


-0 


79 


4 


154 


8 


37 


-G 


20 


-7 


27 





24 


5490 


1480 


140 


60 


4 


870 


22 


00 


-0 


93 


1 


460 


8 


16 


-S 


99 


-7 


56 





26 


£491 


B310 


18 


80 


4 


120 


26 


00 


-2 


19 


4 


186 


8 


22 


-£ 


90 


-7 


74 





00 


5496 


5G20 


13 


90 


4 


370 


26 


00 


-1 


73 


4 


913 


8 


70 


-5 


08 


-7 


35 





27 


5503 


8950 


91 


70 


4 


270 


22 


00 





02 


2 


578 


8 


24 


-G 


05 


-7 


55 





25 


5512 


9800 


141 


70 


6 


140 


20 


00 


-0 


4S 


2 


932 


8 


52 


-4 


05 


-7 


36 





00 


££2£ 


£440 


98 


00 


6 


9£0 


26 


00 


-1 


08 


4 


230 


8 


38 


-£ 


89 


-7 


26 





24 


5520 


8980 


SO 


90 


4 


270 


26 


00 


-2 


02 


4 


473 


8 


22 


-4 


94 


-7 


16 





28 


5529 


IGOO 


58 


20 


3 


930 


26 


00 


-2 


68 


3 


642 


7 


25 


-G 


30 


-7 


59 





24 


5536 


5800 


41 


30 


4 


070 


26 


00 


-3 


73 


2 


832 


7 


26 


-6 


14 


-7 


82 





26 


££43 


9360 


98 


90 


4 


760 


26 


00 


-1 


11 


4 


217 


8 


38 


-£ 


87 


-7 


26 





24 


5546 


5060 


78 


40 


3 


690 


26 


00 


-1 


28 


4 


371 


8 


25 


-5 


40 


-7 


22 





28 


5550 


6090 


24 


60 


2 


890 


72 


00 


-1 


90 





000 





00 





00 





00 





00 


5551 


3330 


7 


70 


2 


840 


41 


00 


-1 


44 





620 





00 





00 





00 





00 



Fig. 4. — Example of the linelist input file (<rootf ilename> . in) for GALA. 



^"W hen one of the 7 constants is 0, WIDTH9 calculates its value adopting some theoretical relations 
2005bl , for more details). 



Castelh 
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6. Output files 

For each star listed in the file list_star, GALA will create three output files. 

(1) <rootf ilenaine> . abu summarizes the derived parameters, the uncertainties and the 
mean abundance for each element available in the input file <rootf ilename> . in. Fig. [5]shows an 
example of this kind of output file. The first thirteen lines are referred to the adopted analysis: 

• the second line is the value of the merit function, that provides an estimate of the goodness 
of the solution; 



• third, fourth and fifth lines concern the correlation among the lines in the A(E1)-EWR , 
A(E1) — Xex and A(E1) — A planes, summarysing the slope, the Jackknife uncertainty in the 
slope, the formal error in the slope, the Spearman rank correlation coefficient and the zero- 
point of the linear fit; 



• the sixth line is the difference between the average abundances of the chosen element derived 
from neutral and single ionized lines, together with its Jackknife uncertainty; 



• the seventh line is a flag (Optimization) about the adopted procedure: 'Optimization 1' 
indicates that the analysis has been performed as specified in autof 1 .param, while 'Opti- 
mization 0' advises that during the analysis the optimization of one or more parameters has 
been switched off. The following three lines specify as log g, Tofj and vt have been optmized. 
The Optimization flag is useful to identify rapidly stars for which the analysis has had some 
problems or failures. 



The final atmospheric parameters are written in the twelfth line, together with the adopted 
model atmosphere (i.e. 'atlasQ', 'marcs' or the name of the used individual model). In case of 
ATLAS9 models, if some atmospheric layers of th e final rnodel atmosphere do not converge 
(according to the convergence criteria proposed by ICastellilll988l ). the number of unconverged 
layers is specified. Also, the label PARAMETERS allows to easily extract the final parameters 
(for instance with grep or awk commands). 



The uncertainties in the parameters are written in the thirteenth line; at the end of the 
line, the paramete rs of the model atmospheres obtained with the procedure described by 
Cavrel et al.l (j2004l ) are reported (only if the option uncer= 1 or 2 is used and Teg has been 



spectroscopically optimized). Note that these uncertainties in the atmospheric parameters are 
calculated by using the Jackknife technique for the optimized parameters: when a parameter 
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is fixed, the uncertainty is that provided by you in list_star, if the latter is 0, the step of 
the grid is taken as representative of the error. 



The following lines provide for each element the mean abundance, the dispersion of the mean, 
the num ber of used liii e s, the net variation of the abundance by using the two models obtained 
with the ICavrel et al.l (|2004l ) method (columns '+Cay' and '-Cay') and the net variation of 
the abundance by varying each time one only parameter (columns '+Teff', '-Teff', '+logg', 
'-logg', '+vt', '-vt'). Additionally, a summary file (named abundance . final) is provided, 
including all the <rootfilename> . abu files for the stars listed in list_star; 
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Fig. 5. — Example of the output file (<rootf ilename> . abu) including information about atmo- 
spheric parameters, uncertainties and average abundances. 



(2) <rootf ilename> . EW includes the main information for each transition included in the 
input file <rootf ilename> . in. Fig. [6] shows a part of this file as reference. The first rows of the 
file explain the meaning of the flag about the line rejection. Briefly: 

flag = 1 the line is used; 

flag = the line is rejected according to the cr-clipping algorithm (see iMucciarelli et al.ll2013l ): 
flag = -1 the line is rejected because outside the range of valid EWR values specified by 'ewmin' 
and 'ewmax' in 'autofl.param'; 

flag = -2 the line is rejected because its error is larger than the boundary value 'error' in autof 1 . param; 
flag = -3 if the abundance calculation performed by WIDTH9 does not converge, producing a 'NaN', 
the value is substituted with 99.999 and the line excluded by the analysis. 
The rest of the flle includes for each transition: 
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• the wavelength (expressed in A); 

• the observed EW (expressed in niA); 

• the logarithm of the oscillator strength; 

• the excitation potential (expressed in eV); 

• the derived abundance; 

• the error of EW, if provided in input (expressed in niA) ; 

• the error in abundance obtained by varying the observed EW of +1(Tew\ 

• the code of the element in the usual Kurucz notation; 

• the flag related to the line rejection; 

• the theoretical EW (expressed in mA) calculated by WIDTH9 through G aussian integra tion 
of the line profile by using the nominal metallicity of the best model (see ICastellill2005bl . for 
details) ; 



• the ne t abundance variation calculated according to the method described by ICayrel et al 



(j2004l ) and keeping fixed the effective temperature at Tgfj + aT^g and T^g — aT^s ; 

the net abundance variations due to the variation of each parameter (ztfiTeff ,ib(Tlog g, ztavt). 
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Fig. 6. — Example of the output file (<rootf ilename> . EW) including information for each indi- 
vidual line. 

(3) <rootf ilename> .ps is the graphical output of GALA where the the main results are 
shown and an example of this output is shown in Fig. [71 The first three panels show the behaviour 
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of the abundances for the chosen element as a function of EWR, Xex and A; the small lower panel 

in the first panel shows the behaviour of the EW error as a function of EWR. The blue lines arc the 
linear fits obtained in each plane. The last panel shows the curve of growth obtained by plotting 
EWR as a function of the theoretical EW^, defined as EWT — loggf - ^Xex- 

Black points are the lines in the ionization stage used to optmize T^s and Vt, while red points the 
lines of the same species but in the other ionization stage; for instance, if you specify 26.00 (Fe I) 
in the 'qiron' parameter, black points in the graphical output will be the neutral iron lines and the 
red points the Fe II lines. Empty points are the lines rejected during the analysis (the reason of 
the rejection is specified in the <rootf ilename> .EW file). 

7. How to manage the grids of model atmospheres 

The automatization of the chemical analysis as implemented by GALA needs the management 
of wide grids of model atmospheres, in order to freely explore the parameter space. Here we briefly 
explain how you can obtain the grids of models. It is worth to recall that you do not need to 
download all the models described in the following: you can only have the grid of models that 
you will use for your analysis. The only caveat is that you cannot call a kind of model in the 
autof 1 .param file if you do not have a directory with all these models. 

If you decide to adopt ATLAS9 models for your analysis, you need to create in the correspond- 
ing directory three sub-directories, named ODF/, MODELS/ and MODELS-new/ (the latter is 
empty and it will be filled with new models progressively using GALA). 

Note that the capability to calculate a given ATLAS9 model atmosphere depends on the 
availability of the ODF computed with the requested chemical composition. The input guess model 
atmosphere has to be as close as possible to the new model in terms of Teg, log g and metallicity in 
order to allow a good convergence, but in principle, the output model is independent from the input 
one. Thus, when a new ATLAS9 model is requested during the analysis process, GALA looks for 
the ODF and Rosseland table corresponding to the requested metallicity, if these tables do exist, 
the closest model with that chemical composition is picked as input guess model (GALA recognizes 
automatically the closest model), otherwise GALA skips the star and advises you with a warning 
message that the selected metallicity is incompatible with the metallicities encoded in GALA. The 
extension (in terms of metallicity, a-enhancement and microturbulent velocity) of each ATLAS9 
grid is written in the source code of GALA, allowing an easy implementation of additional ATLAS9 
grids in the next releases of GALA. 

ATLAS9 (F.Castelli) 

This dataset includes ODFs calculated for 5 values of microturbulent velocities (namely, 0, 1, 2, 
4 and 8 km/s) and two levels of a-enhancement, [a/Fe]=0.0 and +0.4. GALA uses the metal- 
licities between -5.5 and +0.5 dex with a step of 0.5. To easily obtain all the needed files, we 
provide a wget-based procedure. Move the code downl_f x. csh in the directory specified in the file 
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Fig. 7. — Example of the graphical output (<rootf ilename> .ps) of GALA. Left-upper panel 
shows the behaviour of the abundances as a function of the reduced EW, EWR (while the small 
panel shows the behaviour of the EW error as a function of EWR). Right-upper and left-lower panels 
plot the behaviour of the abundances as a function of Xex and of the wavelength, respectively. The 
right-lower panel shows an empirical curve of growth. In all the panels the filled circles are the used 
lines and empty circles are the rejected lines. Black symbols are the lines used in the optimization 
of Tgff and vt (Fe I in this case), while red symbols are the lines in the other ionization stage (Fe II 
in this case). Blue lines are the linear best-fits. 
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GALA/src/Paths.f and type 



. /downljfc . csh 



However, in the following we describe the correct way to obtain and manage the ATLAS9 files. 



1. In the directory specified in the file GALA/src/Paths.f for this grid of models, create three 
subdirectory, named ODF, MODELS and MODELS-new. In the first subdirectory, put ODFs and 
Rosseland tables with [M/H] between -5.5 and +0.5 dex (at step of 0.5) and with both the 
values of a-enhancement. You can download the ODFs from the website of F. Castelli 

http:/ /wwwuser. oat. ts.astro.it/castelli/odfnew. html 

(you need only the big ODFs, do not download the Jit ODFs), as well as the tables with the 
Rosseland mass absorption coefficients 

http: / /wwwuser. oat. ts.astro.it/castelli/kaprossnew. html. 



2. The grids of models are available in the same website 



http : / /wwwuser . oat . ts . astro, it / castelli / grids .html. 

For the metallicities used by GALA, all the models of a given [M/H] are available in a single 
ASCII file (for instance, master file for [M/H]=+0.5 and [a/Fe]= 0.0 is named ap05k2odfnew. dat, 
while for [M/H]=+0.5 and [a/Fe]= +0.4 ap05ak2odfnew.dat.) Save these (14) files in the 
subdirectory MODELS. We recommend to download these files (without renaming them), be- 
cause GALA looks for these files to pick the input guess model (please, avoid downloading 
the grids with peculiar assumptions, as He enhancements and different mixing length values, 
because GALA does not (yet) use these grids to search for the input model). Note that when 
you ask GALA to use a metallicity outside the range -2.5/+0.5, the star will be analysed by 
adopting the closest model. 

The new ATLAS9 models calculated during the runs of GALA will be saved in the subdi- 
rectory MODELS-new, together with a file including the main information about the calcula- 
tion of the model . The new models will be named following the same formalism used by 



Castelli Kuruca (|2004l ) for the individual models; for instance, a model with Tes= 4120 K, 

log g= 1.32, vturb='^ km/s, [M/H]=~0.5 and [a/Fe]= +0.4 wih be saved as am05at4120gl32k2odfnew.dat, 

will its summary file will be named am05at4120gl32k2.1og-at9. 



ATL AS9 (APOGE E) 

Recently, iMeszaros et al.l (j2012l ) computed a new set of ATLAS9 model atmospheres, ODFs and 
Rosseland opacity tables, spanning a huge range of metallicity, a-enhancement and C abundances. 
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In order to avoid the introduction of another input parameter (namely the C abundance, [C/Fe]), 
GALA uses subgrids with fixed C abundances. The models (as well as the ODFs and the Rosseland 
opacity tables) are freely available from the website 

http://www.iac.es/proyecto/ATLAS-APOGEE/ 

Note that if you plan to use only the models with a given [C/Fe] value, you can easily download 
only the desired files or download all the models and ODFs and then delete these not of interest 
(in order save disk space). In following we describe how to obtain all the models but the same 
instructions are valid also to manage ODFs and models for a given [C/Fe] value. 

1. Also for the ATLAS9-AP0GEE grid, create in the corresponding directory three subdirecto- 
ries named ODF, MODELS and MODELS-new. On the webpage, the available data are provided 
in tar.bz2 archive files according to the metallicity. Download the ODFs archive files. 

2. The archive file including all the models (and named alljnod.tar) is available on the main 
page of the ATLAS9-AP0GEE project. Download this file in your subdirectory MODELS/ 
and unpack it. The archive includes for each chemical composition a subdirectory with all 

the individual models and a master file grouping together all the models of that chemical 
composition. Delete all these subdirectories and leave only the master files. These files 
are named according to the chemical composition: for instance, mm05cm03op00 .mod is the 
master file including all the models with [M/H]=-0.5, [C/Fe]=-0.25 and [a/Fe]=+0.0). Do 
not rename neither the models nor the ODFs, otherwise GALA will crash because it will be 
not able to find the needed files to calculate a model atmosphere. 

3. The new ATLAS9 models will be saved in the subdirectory MODELS-new. Also for the 
ATLAS9-AP0GEE models, new models and summary files will be named according to the 

nomenclature used to label the original model of the grids. Also for the models, do not rename 
the files. For instance, the model with Tcff= 4120 K, log g= 1.32, vturb='^ km/s, [M/H]=— 0.5, 
[C/Fe]=-0.25 and [a/Fe]= +0.0 will be named amin05cin03op00t4120gl32v20 .mod. 

MARCS 

MARCS model atmospheres are available at the website 

http://marcs.astro.uu.se/ 

but unfortunately the entire grid of models cannot be downloaded as a whole, because the web- 
interface allows to download a limited (less than an hundred) number of models at a time. If 
you wish to use the MARCS models, he/she needs to patiently download all the models from the 
website and then group them according to their metallicity, microturbulent velocity (as done for 
the ATLAS9 grids) and geometry (plane-parallel or spherical). GALA is designed to manage only 
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the grids with standard chemical composition and microturbulent velocities 1 and 2 km/s, but 
modified versions of the code can be provided on request for other sub-grids of MARCS models. 
Because the grid of models with standard chemical composition naturally varies the a-enhancement 
according to the metallicity, for these models it is not necessary to group them for the value of the 
a-enhancement, at variance with the ATLAS9 models. Each master file needs to have a precise 
name. In Fig. [8] we summarize the name of these master files, as expected by GALA, together with 
the corresponding metallicity and geometry, and the total number of included models. 



marcsjplOOkl 


PP 


[M/H 


= 


+ 1 


00 


vturb 


=1 


km/s 


plane -parallel 


171 


models 


marcs jp 10 0k2 


sp 


[M/H 


= 


+ 1 


00 


vturb 


=2 


km/s 


spherical 


210 


models 


marcs jpOVSkl. pp 


[M/H' 


= 


+ 


75 


vturb 


=1 


km/s 


plane -parallel 


172 


models 


marcs jp075k2 


sp 


[M/H' 


= 


+0 


75 


vturb 


=2 


km/s 


spherical 


219 


models 


marcs jpOSOkl. pp 


[M/H 


= 


+0 


SO 


vturb 


=1 


km/s 


plane -parallel 


174 


models 


marcs jp0S0k2 


sp 


[M/H 


= 


+0 


SO 


vturb 


=2 


km/s 


spherical 


217 


models 


marcs jpOESkl. pp 


[M/H 


= 


+0 


25 


vturb 


=1 


km/s 


plane -parallel 


172 


models 


marcs jp025k2 


sp 


[M/H' 


= 


+0 


25 


vturb 


=2 


km/s 


spherical 


226 


models 


marcs jpOOOkl. pp 


[m/h; 


= 


+0 


00 


vturb 


=1 


km/s 


plane -parallel 


173 


models 


marcs_p000k2 


sp 


[M/H 




. n 

+u 


n n 
00 


vturb 




km/s 


spherical 


Z-iZ 


models 


marcs_m02Skl 


PP 


[M/H' 




-0 


25 


vturb 


= 1 


km/s 


plane -parallel 


173 


models 


marcs_m025k2 


sp 


[M/H 




-0 


25 


vturb 


=2 


km/s 


spherical 


225 


models 


marcs_m050kl 


PP 


[M/H' 




-0 


50 


vturb 


=1 


km/s 


plane -parallel 


173 


models 


marcs_m050k2 


sp 


[M/H' 




-0 


50 


vturb 


=2 


km/s 


spherical 


220 


models 


marcs_m07Skl 


PP 


[M/H 




-0 


75 


vturb 


=1 


km/s 


plane -parallel 


1G8 


models 


marcs_m07Sk2 


sp 


[M/H 




-0 


75 


vturb 


=2 


km/s 


spherical 


211 


models 


marcs_mlOOkl 


PP 


[M/H 




-1 


00 


vturb 


=1 


km/s 


plane -parallel 


168 


models 


marcs_ml00k2 


sp 


[M/H' 




-1 


00 


vturb 


=2 


km/s 


spherical 


210 


models 


marcs_ml50kl 


PP 


[M/H' 




-1 


50 


vturb 


=1 


km/s 


plane -parallel 


162 


models 


marcs_mlS0k2 


sp 


[M/H 




-1 


50 


vturb 


=2 


km/s 


spherical 


202 


models 


marcs_m200kl 


PP 


[M/H 




-2 


00 


vturb 


=1 


km/s 


plane -parallel 


161 


models 


marcs_m200k2 


sp 


[M/H 




-2 


00 


vturb 


=2 


km/s 


spherical 


1S7 


models 


marcs_m250kl 


PP 


[M/H' 




-2 


50 


vturb 


=1 


km/s 


plane -parallel 


165 


models 


marcs_m250k2 


sp 


[M/H' 




-2 


50 


vturb 


=2 


km/s 


spherical 


190 


models 


marcs_m300kl 


PP 


[M/H 




-3 


00 


vturb 


=1 


km/s 


plane -parallel 


136 


models 


marcs_m300k2 


sp 


[M/H 




-3 


00 


vturb 


=2 


km/s 


spherical 


165 


models 


marcs_m400kl 


PP 


[M/H 




-4 


00 


vturb 


=1 


km/s 


plane -parallel 


123 


models 


marcs_m400k2 


sp 


[M/H' 




-4 


00 


vturb 


=2 


km/s 


spherical 


132 


models 


marcs_m£iOOkl 


PP 


[M/H' 




-5 


00 


vturb 


=1 


km/s 


plane -parallel 


108 


models 


^ marcs_m500k2 


sp 


[M/H 




-5 


00 


vturb 


=2 


km/s 


spherical 


117 


models 



Fig. 8. — List of the files with the MARCS models grouped according to their metallicity, micro- 
turbulent velocities and geometry. 
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